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The first excited state belongs to these new levels and displays unusual scaling 
exponents. We conjecture that, while the simple asymmetric exclusion process without 
impurities belongs to the KPZ universality class with dynamical exponent |, our model 
has a scaling exponent | + N— 1. In order to check the conjecture, we solve numerically 
the Bcthc equation with N = 3 and JV = 4 for the totally asymmetric diffusion and 
found the dynamical exponents | and § in these cases. 
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1. Introduction 

The simple asymmetric exclusion model (ASEP) is a stochastic model that describes 
the dynamics of hard-core particles diffusing asymmetrically on the lattice. This model 
became a paradigm in non-equilibrium statistical physics in the same way that the Ising 
model in the equilibrium statistic mechanics. Due to its intrinsic nontrivial many-body 
behavior, the ASEP is used to modeling a wide range of complex systems, like traffic flow 
[1], biopolymerization [2], interface growth [3], etc (see [4] for a review). Remarkably, 
the ASEP in one-dimension is exactly solvable, what enable us to use the Bethe Ansatz 
[5] to obtain spectral information about its evolution operator [6, 7, 8, 9]. The relaxation 
time to the stationary state depends on the system size L and satisfies a scaling relation 
T ~ L z , where z — | is the ASEP dynamical exponent. This dynamic exponent was 
first obtained by the Bethe Ansatz [6, 8, 10] and shows that the ASEP belongs to the 
Kardar-Parizi-Zhang (KPZ) universality class [11]. The scaling property of the model 
can be understood by mapping the ASEP into the particle height interface model, whose 
fluctuations in the continuum limit are governed by the KPZ model [11]. 

On the other hand, the generalization of the simple exclusion problem by including 
more than one kind of particles (N = 1,2,...) has displayed exciting new physics, 
including spontaneous symmetry breaking and phase separation phenomena [12]. The 
introduction of a second class of particle is a useful tool to study the microscopic 
structure of shocks [13], and the case with three distinct classes of particles was first 
considered in [14]. However, the critical phenomena and universal dynamics of these 
one-dimensional driven diffusive systems with several kinds of particles are largely 
unexplored. Another motivation for studying these models stems from the connection 
between interacting stochastic particle dynamics and quantum spin systems. This 
connection follows from the similarity between the master equation describing the 
time fluctuations of these models and the Schrodinger equation in Euclidean time. 
This relationship enables us to identify a quantum Hamiltonian associated for these 
stochastic models. The simplest example is the mapping between the ASEP and 
the exact integrable anisotropic Heisenberg chain, or the so called, XXZ quantum 
chain [4]. Furthermore, iV-state quantum Hamiltonians have played an important role 
in describing strongly correlated electrons in the last decades. Remarkably, in one 
dimension several models in this category are exactly solvable, as for example, the spin- 
1 Sutherland [15] and t-J [16] models, and the spin-| Perk-Schultz model [17], the 
Essler-Korepin-Schoutens model [18], the Hubbard model [19] and the two-parameter 
integrable model introduced in [20]. These quantum models can be related to the 
asymmetric diffusion of two (spin-1) and three (spin-|) kinds of particles [4], respectively. 
In its formulation in terms of particles with two and three global conservation laws, these 
models describe the dynamics of different kinds of particles on the lattice, where the 
total number of particles of each type is conserved separately. 

In order to ensure integrability, all known models of this class satisfy some particle- 
particle exchange symmetries [21, 22]. Recently, we introduced a new class of 3-state 
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model that is integrable despite it do not have particle-particle exchange symmetry [22] . 
In [23] we extend the model [22] and formulate an one-dimensional asymmetric exclusion 
process with one kind of impurities (ASEPI). This model describes the dynamics of two 
types of particles (type 1 and 2) on a lattice of L sites, where each lattice site can be 
occupied by at most one particle. While particles of type 1 can jump to neighbors sites 
if they are empty, like in ASEP, particles of type 2 (called impurities) do not jump to 
empty sites but exchange positions with neighbor particles of type 1. We show that this 
model has a relaxation time longer than the ones for the ASEP, and displays a scaling 
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exponent of z — § [23] (of order L^xL = [23]). We obtained this result by solving the 
Bethe Ansatz equation for the half-filling sector and in the totally asymmetric diffusion 
process [23]. 

In the present work we show how this model can be easily generalized to obtain 
models with relaxation times even larger. We formulate an asymmetrical diffusion model 
oi N = 1,2,3, ... kinds of particles with impurities (N-ASEPI), where particles of kind 1 
can jump to neighboring sites if they are empty and particles of kind a — 2,3, ...,N 
(called impurities) only exchange positions with particles if satisfy a well defined 
dynamics. Different from the ASEPI [23], our generalized model can have more than 
one particle on each site (multiple site occupation). Although our model can be solved 
by the coordinate Bethe ansatz, we are going to formulate a new matrix product ansatz 
(MPA) [24, 21] due its simplicity and unifying implementation for arbitrary systems. 
This new MPA introduced in [24, 21] can be seen as a matrix product formulation of 
the coordinate Bethe Ansatz and it is suited to describe all eigenstates of integrable 
models. We solve this model with periodic boundary condition through the MPA and 
we analyze the spectral gap for some special cases. Our N-ASEPI model displays the full 
spectrum of the ASEP [6] plus new levels. The first excited state belongs to these new 
levels and displays unusual scaling exponents. Although the ASEP belongs to the KPZ 
universality class, characterized by the dynamical exponent z — § [11], we conjecture 
that our model displays a scaling exponent | + N — 1, where N — 1 is the number of 
kinds of impurities. In order to check our conjecture, we solve numerically the Bethe 
equation with N = 3 and N = 4 for the totally asymmetric diffusion and found that 
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the gap for the first excited state scales with L~z and L~a in these cases. Furthermore, 
we also generalize the model [23] to include quantum spin chain and solve the Bethe 
Ansatz equation for symmetric and asymmetric diffusion. 

Our paper is organized as follows. In section 2 we generalize the model [23] to 
include quantum spin chain and solve the Bethe Ansatz equation for the symmetric 
and asymmetric diffusion. The generalization for several kinds of impurities is done in 
section 3. Finally, our conclusions are presented in Section 4. 

2. The asymmetric exclusion model with one kind of impurities 

Recently, we propose an exactly solvable asymmetric exclusion process with impurities 
[23] (ASEPI) and found its dynamic exponent z — §. The exponent z in [23] was 
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obtained, from the spectral gap of the model, for the totally asymmetric exclusion 
process (TASEPI) and at half-filling. It is important to notice that although the ASEP 
without impurities belongs to the KPZ universality class [11] (dynamic exponent |), 
our new model displays an unusual scaling exponent |. In this section we extend our 
previous analysis [23] and obtain the spectral gap for the symmetric and asymmetric 
exclusions process. Furthermore, we generalize both the models [22] and [23] in order to 
include quantum spin chains, and we found analytically the spectral gap of the quantum 
model in the special case where we have free fermions. 

The model in [23] describes the dynamics of two kinds of particles (type 1 and 
2) on an one-dimensional lattice of L sites, where each lattice site can be occupied by 
at most one particle. Furthermore, the total numbers ni, n 2 of particles of each type 
is conserved. In this model if the neighbor sites are empty, particles of type 1 can 
jump to the right or to the left with rate Tq5 and T^q, respectively. Particles of type 
2 (impurities) do not jump to neighboring sites if they are empty, but can exchange 
positions with neighbor particles of type 1 with rates \ and Y\ \ if the the particle 1 
is on the left or on the right, respectively. To describe the occupancy of a given site i 
(i = 1,2, ...,L), we attach on it a variable taking values cui = 0, 1,2. If «j = 0, the 
site is vacant. If ai = 1, 2, we have on the site a particle of kind 1 or 2, respectively. The 
allowed configurations can be denoted by the set {a} = {a±, a 2 , «l} of L integers 
on = 0,1,2. The master equation for the probability distribution at a given time t, 
P({a},t), can be written in general as 

dn{ ^ t] = nW} {*})P(W,t) - r(W W})P({a},t) (i) 

where r({a} — > {a'}) is the transition rate where the configuration {a} changes to {a'}- 
The master equation (1) can be written as a Schrodinger equation in Euclidean 
time (see [4] for general applications for two-body processes) 

where we represent a configuration oii on site i by the vector [a^j, and we interpret 
\P) = P({a} : t)\cti) <g> [0:2) <8> • • • <E> \chl) as the associated wave function. In order 
to generalize our model [23] and to include quantum chains solutions, the general 
Hamiltonian we consider on a ring of perimeter L is given by: 

3=1 \ a#3=l 

+ E K^E^y (3) 

a,P=Q J 

with = E^ due to the periodic boundary condition, and where E^ ,lB (a,/3 = 

0, 1, 2) are the 3x3 Weyl matrix acting on site k with i,j elements [E^j = 5i t i5 m j, 

and T^™ are the coupling constants. The last sum in (3) accounts for the static 
interactions while the first and second sums are the kinetic terms representing the motion 
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and interchange of particles, respectively. The U(l) <S> U(l) symmetry supplemented 
by the periodic boundary condition of (3) imply that the total number of particles 
ni, n 2 = 0, 1, 2..., L (with n\ + n 2 < L) on class 1 and 2 as well the momentum P = ^ 
(I = 0, 1,...,L — 1) are good quantum numbers. Furthermore, the Hamiltonian (3) 
also preserves the numbers of vacant sites between the impurities. This conservation 
plays a fundamental role in the spectrum properties of the model [23]. As we shall 
show, for the stochastic model the Bethe equation do not depends on the number of 
impurities (n 2 ^ 0). Consequently, the roots of Bethe equation and the eigenvalues 
of the Hamiltonian are independent of n 2 ^ (but the wave function depends on n 2 ). 
This huge spectrum degeneracy follows directly from the conservation of the numbers 
of vacant sites between the impurities by the Hamiltonian (3). Let us explain with 
the following example. Suppose we start with a given configuration 0120220 with one 
particle (1), 3 impurities (2) and 3 vacant sites (0). We can make a surjective map 
between all possible configuration of these particles to all possible configurations of a 
new chain with just impurities and vacant sites. For example 0120220 =>- 020220. On 
this new chain, we are looking only for the effective movement of impurities on the 
chain. For simplicity, let us consider the totally asymmetric model (TASEP) where 
T\ \ = 1 and T^q = 0. When the particle jumps over the impurities, nothing changes 
in the effective chain since we also have 0210220 =>• 020220, then 0201220 =>- 020220, 
then 0202120 020220, then 0202210 020220, then 0202201 020220, and 
finally a change in the mapped configuration 1202200 ==>- 202200. On other words, 
the impurities move on the mapped chain as they are just one "object" due to the 
conservation of vacant sites between impurities. Moreover, this "object" only moves 
when the particle complete a turn over the chain. As a consequence, the time for the 
particle to complete one revolution is the time scale for the movements of the this 
"object". For an arbitrary number of particles in a chain of length L, the time for the 
particles to complete one revolution is of order (L 2 in the symmetrical diffusion). 
As the "object" formed by the impurities need to move of order L times to span all 
possible configurations, it will takes a time of order L? x L = La (L 3 in the symmetrical 
diffusion) to reach the stationary state. 

2.1. The exact solution of the model 

We want to formulate a matrix product ansatz for the eigenvectors \^ ni ,n 2 ,p) of the 
eigenvalue equation 

belonging to the eigensector labeled by (ni,n 2 ,P). These eigenvectors are given by 

\^n u n 2 ,p) = J2J2f(xi,ai, . . . ;x n ,a n )\xi,ai\ . . .;x n ,a n ), (5) 

where the kets \xi, ai, . . . ; s n , a n ) = (\0)^>) Xl ~ 1 \ai) <S> (\0)^>) X2 ^ Xl ^ 1 \a 2 ) <%>•••<%> 
(lO)®)^"^ 1-1 \a n ) ® (\0)®) L ~ Xn denote the configurations with particles of type a { 
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(aii — 1, 2) located at the positions Xi (xi — 1, . . . , L), and the total number of particles 
is n = U\ + ri2- The summation {a} = {cti, . . . , a n } extends over all the permutations 
of n integers numbers {1, 2} in which n\ terms have value 1 and n2 terms the value 2, 
while the summation {x} = {x\, . . . ,x n } extends, for each permutation {a}, into the 
set of the non-decreasing integers satisfying > Xi + 1. 

The MPA [24] is constructed by making a one-to-one correspondence between the 
configurations of particles and a product of matrices 

f{x u ai ;...;x n , a n ) <^ e^A^E^'^A^ ■■■ (6) 

£] x n ~ x n - 1 - 1 A («n ) £jL - X n 

where for this map we can choose any operation on the matrix products that give 
a non-zero scalar. In the original formulation of the MPA with periodic boundary 
conditions [24] the trace operation was chosen to produce this scalar. The matrices 

are associated to the particles of type a = 1,2, respectively, and the matrix E 
is associated to the vacant sites. Actually E and A^ are abstract operators with an 
associative product. A well defined eigenfunction is obtained, apart from a normalization 
factor, if all the amplitudes are related uniquely, due to the algebraic relations (to be 
fixed) among the matrices E and A^K Equivalently, the correspondence (6) implies 
that, in the subset of words (products of matrices) of the algebra containing n matrices 
AW and L — n matrices E there exists only a single independent word ("normalization 
constant"). The relation between any two words is a c number that gives the ratio 
between the corresponding amplitudes in (5). 

As the Hamiltonian (3) commutes with the momentum operator due to the periodic 
boundary condition, the amplitudes f(xi, a±; . . . ; x n , a n ) should satisfy the following 
relations: 



where 



f(x 1 , a 1 ;...',x n ,a n ) = e lP 'f{x 1 + 1, a^, . . . ; x n + 1, a n ), (7) 

P=^, / = 0,1,...,L-1. (8) 

Let us consider initially the simpler cases where n — 1 and n — 2. 

n = 1. We have distinct equations depending on the type a — 1, 2 of the particle. 
The eigenvalue equation (4) give us 

£ (l) E x-l A (l) E L-x = r 10 E x-2 A (l) E L-x+l + pO l E x A (l) E L-x-l 

+ (t° 1 + t\ )e x ~ 1 a^e l - x , (9) 

if the particle is of type 1 and 

£ (2) E x-l A (2) E L-x = ( r 2 + T 2 0) ^-1^(2)^-^ (1Q) 

if the particle is of type 2. In these last two equations = e 1 ' and = e 0,1 are the 
eigenvalues, and we choose ° = without loss of generality. A convenient solution is 
obtained by introducing the spectral parameter dependent matrices 

AW = EA£> (a = 1,2), (11) 
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with complex k parameter, that satisfy the commutation relation with the matrix E 

EA^ = e lk A ( k a) E (a = 1,2). (12) 
Inserting (11) and (12) into (9) and (10) we obtain 

e^Cfe) = rS ?e-* + r? ie* + rg i + ri g, 

= rg 2 + r28. (is) 

The up to now free spectral parameter k is fixed by imposing the boundary condition. 
This will be done only for general n. 

n = 2. For two particles of types a\ and ct2 (011,0:2 = 1,2) on the lattice we have 
two kinds of relations coming from the eigenvalue equation. The configurations where 
the particles are at positions (xi, X2) with X2 > x\ + 1 give us the generalization of (9) 

£ ni, 712 jjjxi—l j^(ai) jjjxz—xi —1 j^az) jjjL—X2 _ 

pai jj*x\—2 J^( a i) £X2—xi j^(a 2 ) j^L-x 2 _|_ pO a\ jjjxi j^i^) jj<X2-xi-2 j^(a 2 ) £jL~X2 
_|_ p"2 £JX!-1 ^(«i) j^x 2 ~x 1 -2 j^(a 2 ) j^L-x 2 +l 
_|_ pO a2^ia;i-l^(Q:i) JjJX2—Xl j^(oi2) ■£jjL—X2 — 'L 

_|_ ^pO ai _|_ pai _|_ pO 0L2 _|_ pQ2 0~j j^xi-1 J^( a i) j^x 2 -xi~l j^(a 2 ) j^L—x 2 (14) 

and the configurations where the particles are at the colliding positions (xi = x, 
X2 = x + 1) give us 

£.ni,n,2 £jx— 1 j^(ai) j^(oi2) j^jL— x— 1 po-i jjjx— 2 j^(a-i) j^(ct2) j^L—x—l 

_|_ pO 0-2 j^x~ ^ J^{oti) j^(a 2 ) j^L— x— 2 _|_ p«2 oi\ j^x—1 j^(a 2 ) j^(cti) j^L—x—1 

+ (rg Z +Kl°o + K\ Z) e^a^a^e^- 1 , (15) 

where we introduced Tq^ = T^g = 0, and °i — if oi\ — a^. The Hamiltonian (3) 
do not have a standard solution as in [24] where each of the matrices A^ (a = 1, 2) 
are composed by two spectral parameter matrices, with the same value of the spectral 
parameters ki,k 2 (case a in [24]). In order to obtain a solution for (14)-(15) we now 
need to consider the A^ as composed by n a spectral parameter dependent matrices 
A^ a) belonging to two distinct sets of spectral parameters [22, 23], i. e., 



AM=J:eaM with EA^^e^A^E, (a™ )* = <), (16) 

for a = 1,2 and n\ + n2 = n. These last relations when inserted in (14) give us the 
energy in terms of the spectral parameters k^ (a — 1, 2) 

ni ri2 

e n ^ = J2e {l \kf) + J2e {2) (kf), (17) 

3=1 3=1 

where e^ a \k) is given by (13). 
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Let us consider now (15) in the case where the particles are of the same type. For 
two particles, when ot\ = a 2 = 1, the equations (16), (17) and (15) implies that the 
matrices {A^} should obey the Zamolodchikov algebra [25] 

A W -QiifuW t(i)u(i) aW (aj.i\ ( aW 



^A^^Sl^^knA^A^ (j^l), IA"\ =0, (18) 

3 I I 3 \ 3 ' 

where j, / = 1, rii, and the algebraic constants SlKkj 1 ^, k^) are given by: 



pi , pO lj(k) >+k\ ') _ (r i l _ pi o _ F 1\ ' 
qiifi.W kM)\ — l Z 1 I 1 1 1 1 1 1 lJ e - /-i q\ 

111 j ' 1 J " r io + r oV( fe « +fe «)_ rr ii_ r io_ r oi le ^- Uyj 

1 1 + 1 1 e J I 1 11 i 10 x l^ 1 e ' 

For two impurities (cci = a 2 = 2) at "colliding" positions, the eigenvalue equation does 
not fix a commutation relation among the matrices A {2) since (15) is automatically 

satisfied in this case. On the other hand, the sum 

E A % E * A % * 0. (20) 

where the number of vacant sites between the impurities d = y — x is a conserved 
charge of the Hamiltonian (3), should be different from zero or the MPA will produces 
an eigenfunction with null norm. Moreover, the algebraic expression in (16) assures 
that any matrix product defining our ansatz (6) can be expressed in terms of two single 
matrix products A^ 2) A^ 2) E L and A^ 2) A^ 2) E L . Using (16) we have, from the periodic 

boundary condition, 

A%A%E<- = e^?\-^A%A$E\ (21) 

3 I 3 I 

(2) (2) 

To satisfy this equation we should have k 2 = —k\ + 2nj/L (j = 0, 1, L — 1). Con- 
sequently, the most general commutation relation A^ 2) A^ 2) = Sff (/c,- 2 ' ) , k^^A^A^ 

k-^ ^2 ^2 ^1 

(2) (2) (2) (2) (2) (2) 

among the matrices A (2) and A (2) can be reduced to A (2) A (2) = A (2) A (2) 

k-^ /c 2 ^1 ^2 ^2 ^1 

(5*1 K^j 2 '*, fcj 2 ^) = 1) by an appropriate change of variable in the spectral parameter 
kf\ By choosing S^ikf 1 ikf^) = 1 and imposing that the sum (20) is not zero, we 
obtain (j, I, v — 1, n 2 ) 

*p ^ *» + 5^±i) (m = , !,...), (22) 

where {d v } is the set of all numbers of vacant sites between the impurities. 

Let us consider now the case where the particles are of distinct kinds. From (13), 
(16) and (17), equation (15) give us two independent relations: 



-pa 2 I -pO ai tffj^'+fc' 2 ') / -pcti a 2 pai pO a 2 \ 

1 a 2 "T 1 «i e ^ Qi a 2 "1 1 a 2 / 



fc(«2) 



/1 fc( a i) /1 fc( a 2) 



-rsSe* (aa) ^)4?-!) = ° («i^« 2 = l,2). (23) 

This two relations need to be identically satisfied, and since at this level we want to 
keep k^ and k^ as free complex parameters, (23) imply special choices of the coupling 
constants T™ n [22, 23]: 

r° o = r 2 ) ° — 0; 2^1 i = 1^1 ?> *i2 — ^21 = ^22 = 0, (24) 
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where t aia2 = F%\ % - g - rg^ (ai,a 2 = 1,2). We also obtain the structural 
constants: 

The integrability conditions (24) generalizes the result obtained for the stochastic 
process with one kind of impurities [23] to include quantum chains [22]. Let us consider 
now the case of general n. 

General n. We now consider the case of arbitrary numbers ni, n 2 of particles of 
type 1 and 2. The eigenvalue equation gives us generalizations of (14) and (15). To 
solve these equations we identify the matrices as composed by n a spectral dependent 
matrices (16). The configurations where x i+ i > Xi + 1 give us the energy (17). The 
amplitudes in (15) where a pair of particles of types a,\ and 0:2 are located at the closest 
positions give us the algebraic relations 

j K i K i K j 

where the algebraic structure constants are the diagonal ^-matrix defined by (19), (25), 

S%l(kf\kf :) ) = 1, with coupling constants (24). 

In order to complete our solutions we should fix the spectral parameters 
. . . , fcW and kf\ . . . , fcg). The algebraic expression in (16) assures that any matrix 

product defining our ansatz (6) can be expressed in terms of the matrix product 
■ ■ ■ A^A^ 2) ■ ■ ■ A^ 2) E L . From the periodic boundary condition we obtain: 

&i k n ^ k>%2 



V 1 1 0/ i=i 

e^ )(L - Bl) = fSfV 1 (J = l,-.,n 2 ), (27) 



where the spectral parameters {£;] } should satisfy the restrictions (22). We can write 
the Bethe equation (27) in a more convenient way. From the second expression on (27) 
we have 

n ~L n 2 

e .E Sl f<w, = (nj)™- ^ e .E;:, *f = (gj) ^ (28) 

with m — 0, 1, L — n\ — 1. By inserting (28) and using (19) in the first equation in 
(27) we obtain 



(1) ni 



T^ + T%e l(k y' +k ^-2Ae ik y 



where j = 1, ...,ni, 2A = r[j ° + \ — Tj ° — T[j J, and the phase factor 4> m ,n 2 ( m ) is 
defined by: 




0n 1 ,„ 2 (m)= e -^-i m (m = 0,l,...,L-m-l). (30) 
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The Bethe equation (29) generalizes our results [22, 23]. Furthermore, it is important 
to notice that (30) differs from the ones related to asymmetric XXZ chain [6] by the 
phase factor <f> ni ,n 2 { m ) (30). As we shall see, this phase factor will play a fundamental 
role in the spectral properties of the model. 

Finally, the eigenstate momentum is given by inserting the ansatz (16) into the 
relation (7): 

P = T,kf ) + T,kf = -^ (/ = 0,1,...,L-1). (31) 

3=1 3=1 L 

The Bethe equation (30) plus the momentum equation (31) completely fix the spectral 
parameters {kf^} and {kf^} and the eigenvalues (17). Let us consider some special 
cases: 



2.2. Stochastic Model 

For stochastic models we should have f = — ^ (a, (3 = 0, 1, 2). Let set, without loss 
of generality, ° + T° J = 1, T\ \ = Y\ J and A = |. In this case our model describes an 
asymmetric exclusion process with impurities [23]. The Bethe equation (29) with the 
phase factor (30) reduces now to 



" ( ' 1 1 irj? + r?s e <'«.'">-e«=: 1 " ( ' 

where j = 1, n\, and m — 0, 1, L — n\ — 1. 

In our previous work [23] we consider the TASEPI (when Y\ \ = 1 and J = 0, or 
rj 5 = and r°Q = 1). In this case we solve (32) numerically up to L = 1024 in the half- 
filling sector ni = L/2 and we obtain the scaling exponent z — | for the TASEPI. Now, 
we generalize our previous results by solving the Bethe equation (27) for the asymmetric 
ASEPI and symmetric SEPI exclusion process. We also checked the eigenvalues obtained 
from exact diagonalization of the Hamiltonian with the Bethe Ansatz solution for a small 
chain with L = 6, n\ = 2, n 2 = 1 and To? = 0.75 (see Appendix A). The eigenvalue 
with the largest real part is e ni ' n2 = corresponding to the stationary state (it is 
provided by choosing m = and P = in the Bethe equation (32) giving us the n\ 
fugacities e l i =1)- Others eigenvalues contribute to the relaxation behavior to the 
stationary state. In special, the eigenvalue with the second largest real part determines 
the relaxation time and the dynamical exponent z. This eigenvalue is obtained from 
the Bethe equation (32) by choosing m = 1 and P = ^ (see Appendix A for a detailed 
discussion for a small chain). In Table (1) we show the dynamical exponent z versus Y\ \ 
obtained from the numerical solution of the Bethe equation (27) for several values of L 
in the half-filling sector n\ = L/2. The errors displayed are computed from the linear 
regression for the logarithm of the real part of the energy gap versus the logarithm of 
L (see Figure (1) for the cases where T^® = 0.7 and Tq® = 0.5). In the ASEPI, we 
consider the cases where Tj? = 0.9,0.8,0.7,0.6 with L = 20,40,80,160,200,300,400, 
and we found that the energy gap has a leading behavior of KPZ and a sub-leading 
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Ln(L) 



Figure 1. In this figure we display the logarithm of the real part of the energy gap 
versus the logarithm of L for ° = 0.7 and rj ° = 0.5. The dynamical exponents z 
are 2.50(2) and 3.01(2), respectively. 

Table 1. Dynamical exponent z versus ° in the half-filling sector rii = L/2 





rj? = o.9 


rj° = o.8 


Tl ? = 0.7 


rj? = o.e 


? = o.5 


z 


2.52(2) 


2.51(2) 


2.50(1) 


2.52(2) 


3.01(2) 



term L~z related to the super-diffusion of particles 1 and sub-diffusion of particles 2, 
respectively. For m = (and also for n2 = 0) the sub-leading term vanishes and we 
recover the spectrum of the ASEP without impurities [6]. For m = 1 the spectral gap 

5 3 

scales with L~a instead of due to the sub-leading term. On the other hand, we 
obtain for the SEPI (rj? = T ^ = 0.5 with L = 20,32,40,52), a energy gap with a 
leading behavior of L~ 2 and a sub-leading term L~ 3 related to the the normal diffusion 
of particles 1 and sub-diffusion of particles 2, respectively. 

Finally, as in [23], it is important to notice that the dynamics of particles 1 is not 
affected by the impurities. This explain why the model displays the full spectrum of the 
ASEP [6]. The dynamics of the impurities is totally dependent on the particles, since 
the impurities only move when particles change position with them. Consequently, the 
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time to vanish the fluctuations on the densities of particles acts as a time scale for the 
diffusion of the impurities, resulting in a relaxation time greater than the one for the 
standard ASEP and SEP (reflected in the x L^ 1 = gap for the ASEPI and 
L~ 2 x L- 1 = L- 3 gap for the SEPI). 

Eigenstates: The stationary state is the eigenstate associated to the eigenvalue 
£ m,n 2 _ q ^ j g m j- ne sec tor with m = and P = 0). In this case, the Bethe equation 

, (a) 

(32) has the unique solution e ' =1 for all a = 1, 2 and j = 1, n a . As a consequence, 
the S'-matrix reduces to the identity S°^^(kj ai \ki a2 ^) = 1, and all amplitudes in the 
eigenfunction (5) becomes equals to a normalization constant f(xi, a±; . . . ; x n , a n ) = fo, 
since we have from (6) and (16): 

E X1 A { $E X2 ~ X1 A { $ ■ ■ ■ E^-^A^E 1 -^ (33) 

- A {1) ■■■ A {1) A {2) ■■■ A {2) E L 

The eigenfunction |^ ) corresponding to the stationary state is given by a simple 
combination of all possible configurations of particles, where each particle configuration 
has the same weight given by the normalization constant f . We have from (5): 

l^o) = foJ2J2 kijai; • • • ;x n ,a n ). (34) 

{a} {x} 

Consequently, at the stationary state, all configurations that satisfy the hard-core 
constraints imposed by the definition of the model occur with equal probabilities given 
by f . Furthermore, in this case each site is occupied by a particle a with probability 
p a = n a /L. 

Finally, the MPA (6) enable us to write all eigenstates of (3) in a matrix product 
form. For a given solution and kf\ the matrices E and A^ a) have the following 

finite-dimensional representation: 

E = ® ( l ^> I I n | , (35) 

A $> = ®( o r -> ) ^ /2 ^( i o ) ® /2 ' 

3 1=1 \ U fYi e / 1=1 \ 1 U / 1=3+1 

ni+ri2 

h 

3=3+1 

where n — ri\ +n2, 12 is the 2x2 identity matrix, and the dimension of the representation 
is 2 n . The matrix product form of the eigenstates are given by inserting the matrices 
(16) defining the MPA (6) with the spectral parameters (29) into equation (5). As 
showed in [27], our MPA generalizes the steady-state Matrix Product introduced by 
Derrida et al [28]. For the stochastic model the stationary state is obtained by choosing 
k^p = and k^ = in (35). However, a relation between our matrix product form for 
the steady-state and the standard Matrix Product form for open boundary condition is 
not trivial [27, 28]. 
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2.3. Quantum model 

We consider here only the simplest case of free fermions, when A = and Y\\ = 
i — T?o = 1, in this particular case, the Bethe equation (29) reduces to roots of unit. 
The simplicity of Bethe equation enable us to calculate analytically the eigenvalues by 
following [6]. As in the stochastic model, the ground state is obtained by choosing m = 
and the first excited state is given by m — 1. Due to the L _1 term in the phase factor, 
the energy gap 



with a = 1(2) for ri\ even (odd), scales with L~ 3 instead of L" 1 for the case of only 
one kind of particle. The boundary condition plays a fundamental role in the scaling 
behavior of the model. In the quantum sector our model is related to the strong regime 
of the t-U Hubbard model introduced in [26] and solved with diagonal open boundary 
condition. Although the open chain t-U Hubbard model has a scaling gap of L _1 , since 
its spectrum coincides with the spectrum of the anisotropic XXZ model [26] at A = 0, 
our model with periodic boundary condition displays a scaling gap of L~ 3 . 

3. The asymmetric exclusion model with A — 1 kinds of impurities 

We generalize the model discussed in the previous section by adding more types of 
impurities. Although we can formulate a general model including both stochastic 
process and quantum spin chains, we will consider for simplicity only stochastic process. 
The model introduced describes the dynamics of N types of particles on an one- 
dimensional lattice of L sites, where the total number ni, n^, of particles of 
each type is conserved. Different from the case with just one kind of impurity [23], 
discussed in the previous section, in our generalized model we can have more than one 
particle on each site (multiple site occupation). In order to describe the occupancy 
of a given site i (i — 1,2, ...,L) we attach on site i a set {a}i = {a±, a n }, where 
ctj = 1,...,A (j = l,...,n) denotes a particle of kind otj. If {a}i = 0, the site is 
vacant. If {a}i = {a±, ...,a n }, we have on the site n particles of kinds a±, ct2,..., 
a n with (Xj + i > ctj + 1 (j = l,...,n — 1). The allowed configurations, denoted by 
{a} = {{a} i, {a}2, {«}l} are those satisfying the hard-core constraints imposed by 
the condition o^ +1 > a, • + 1 for particles on the same site (see Figure (2) for an example 
of an allowed configuration). 

The master equation for the probability distribution at a given time t, P({a},t), 
can be written in general as 



where r({o;} — > {a'}) is the transition rate where a configuration {a} changes to {a'}. 
In the model proposed there are only diffusion processes. As in [23], and in the last 
section, if the neighbor sites are empty, particles of type a — 1 can jump to the right or to 




(36) 



dP({a},t) 



T({a} K})P({a}, t) + T({a'} {a})P({a'} : t) (37) 



dt 
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Figure 2. In this figure we have an allowed particles configuration for L = 9, N = 5 
and n = n\ + + n 3 + m + n 5 = 9. The configuration satisfy the hard-core constraints 
imposed by the condition a^+i > (Xj + 1. 



the left with rate Tq 5 and J, respectively. Particles of type a = 2, N (impurities) do 
not jump to the neighbor sites if they are empty. The only allowed motions for impurities 
are those in which 21 (I = 1, 2, ...) particles exchange positions (when ai + i > 21 + 1) 

{1, 3, 21 - 1, a l+1 , ...h{2, 4, 2/, A+i, -h+i (38) 

{2, 2/, a, +1 , ...},,{1, .., 2Z - 1, A+i, - W 

with transition rate Tj^, and 

{2, 21, /3 l+1 , ...}<{1, .., 2Z - 1, a m , (39) 

=S> {1, 3, 2Z - 1, A+i, ...}i{2, 4, 2Z, a m , ...} m 

with transition rate 1^ J, respectively, or those in which 2Z+1 particles exchange positions 
(when > 21 + 2) 

{1, 3, 2/ + 1, a i+2 , -}i{2, 4, 2Z, A+i, - W (40) 

{2, 2Z, a /+2 , 2Z + 1, A+i, 

with transition rates Tq5, and 

{2, 21, ...},{1, 2Z + 1, a l+2 , (41) 

=S> {1, 3, 21 + 1, A+i, ...}i{2, 4, 2Z, a l+2 , ...} m , 

with transition rates J, respectively. It is important to notice that all allowed motions 
are those in which we have a sequence of particles 1,3, ...,2Z ± 1 on site i (i + 1) and 
a sequence of particles 2,4, ...,2Z on site i + 1 (i), respectively (see Figure (3) for an 
example of an allowed motion of particles). 

The stochastic Hamiltonian associated to the master equation (37) on a one- 
dimensional lattice of L sites and periodic boundary condition is given by 

H=t [rSSC^ifi - eI'e^ + r°l(Eh°Ek\ - i^Xii) (42) 

1=1 

, r a (3 / jp{P'}i,{a}i TpW}i+i,{P}i+i T7'{a}i,{a}i rp{P}i+i,{P}i+i 

+ 2^ l p'a'\- tj i ~ ^t+l 

{a}i,{/3} i+1 
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Figure 3. In this figure we have an allowed motion of four particles with transition 



rates rJJ. 



where the summations {a}i and extends over all possible configurations satisfying 

the hard-core constraints imposed by the condition atj+i > atj + 1 and fa + % > fa + 1 
for particles on the same site, and E\ a} " {Ph (with = E[ ahmi due to the 

periodic boundary condition) are operators that annihilate the configuration {fai and 
create a configuration {a}i on site i. The coupling constants T%^ a , in (42) are equal to 
rj° or r?J if {/3} i+1 -> {/?'}; satisfy the relations (38) and (39), or (40) 

and (41), respectively, and are zero otherwise. 



3.1. The exact solution for the asymmetric exclusion model with N — 1 types of 
impurities 

As in the case of one kind of impurity discussed in the previous section, we want to 
formulate a matrix product ansatz for the eigenvectors | 1 ^ m ,na I .»,njv J i') °f the eigenvalue 
equation 

H\y ni , n2 ,..., nN , P ) = e ni ' m '-' nN \^ ni ,n 2 ,...,n N! p) (43) 

belonging to the eigensector labeled by (ni, n 2 , ...,n N , P). The MPA we propose asserts 
that the amplitudes corresponding to the configurations where there are no multiple 
occupancy are given by 

/(si,oi;...;s w Q! ft )^ (44) 

while if there exists a multiple occupancy with m particles fa, fa, f3 m (fa+i > fa + 1) 
at Xj we have 

/(...; xj, fa; xj, fa; ...; x j} /3 m ; . . .) (45) 

where the matrices are associated to the particles of type a (a = 1,2,3), the 
matrices — A^E -1 A ^ E~ l ■ ■ ■ A^ m ^ are associated to a multiple occupation 
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of m particles /3 2 , An at same site, and the matrix E is associated to vacant sites. 
Furthermore, as the Hamiltonian (42) commutes with the momentum operator due to 
the periodic boundary condition, the amplitudes f(x±, ai, . . . ; x n , a n ) should satisfy the 
following relations: 

f(x u ai, . . . ; x n , a n ) = e~ lP f(x 1 + 1, «i; . . . ; x n + 1, a n ), (46) 

where 

P=—, z = ,l,..,L-l. (47) 

The eigenvalue equation (43) give us two kinds of relations for the amplitudes (44) 
and (45). The first kind is related to those amplitudes without multiple occupancy, and 
the second type is related to those amplitudes with multiple occupancy. Let us consider 
separated each case: 

Without multiple occupancy: In this case, we have from the eigenvalue equation 
(43) relations for amplitudes without collisions (particles of kind 1 have only empty 
neighboring sites) and with collisions. For configuration without collision the amplitudes 
should satisfy the following constraints: 

e n ^-' nN f(x 1 ,a 1 ;...;x n ,a n ) (48) 

n 

= E 'i'o,,".A-:-'-- - «*; •••) + r°f6/(-; *i + I, ",: •••)' 

1=1 

- n 1 f(x 1 ,a l ; ...;x n ,a n ), 

where we introduced r° q = Tq ° = for a ^ 1. In order to obtain a solution for 
(48) we need to generalize (16) for iV kinds of particles. We consider the matrices A^ 



to N distinct sets of spectral parameters, i. e., 



(a = 1, N) as composed by n a spectral parameter dependent matrices A ( ^ belonging 



A^=J2EA^ with EA^^e^A^E, (^) =0, (49) 

for a = 1, 2, JV and n± + n 2 + • • • + un — n. These expressions when inserted into 
relations without collisions (48) give us the energy in terms of the spectral parameters 

e ni ' n2 '-' nN =J^e{kf y ), (50) 



where 



e(jfe) = Tl<ie- lk + T<ile lk - I. (51) 

It is important to notice that, like in the previous section, the eigenvalues of the model 
depend only on the spectral parameters of particles of kind 1. On the other hand, for 
amplitudes without multiple occupancy and particles of kind a and (3 at the colliding 
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positions (xj+i — Xj + 1), the eigenvalue equation (43) give us the generalizations of 
(26) with a, — 1, 2, where the S'-matrix elements are given by 

cII/lW i 01+ i 10 e /rr,X 

1 o i ~r 1 i o e J — e' 

-2 1/. (2) , _ 1 _ p ifc (2) 



^K^ 1 ),^ 2 )) 

We do not consider here the cases where we have one particle of kind 1 and other particle 
of kind greater than 2 since the eigenvalue equation in this case will relates amplitudes 
with multiple occupancy. For a, (3 > 2 the relations coming from the eigenvalue equation 
are identically satisfied. In this case, as in the previous section, we can choose without 
loss of generality 

S::(kf,ki a) ) = l (a>2). (53) 

Furthermore, in order to have amplitudes with non null norm, the spectral parameters 
{kj*^} (a > 2) should satisfy the constraints 



t ] * kf + ni2 ™ + h (m = 0, 1, ...) (a > 2), (54) 



where {d^} is the set of all numbers of vacant sites between particles of type a > 2. 

With multiple occupancy: Let us consider first the relations coming from the 
eigenvalue equation (43) where we have only one site with multiple occupancy of m 
particles f3±, j3 2 , fa (A+i > A + 1) at position Xj = x. We have from the eigenvalue 
equation (43) the following relations: 

e n ^- n »f(..-x,p i; x,(3 2 ;..-x,p m ;...) (55) 
= E I'ii ..".A --^ — •••) • rl! "V,/(. • l.n,:...) 

+ r^J(..-x-lJ 1 -,x,(3 2 -...,x,f3 m -,...) 

+ r Sfo/(-; X > fa •••> X > Pm] X + l, Pi, ...) 

-n 1 f(...;x,/3 1 ;x,/3 2 ; x, /3 m ; ...), 
for empty neighboring sites, and 

£ n u n 2 ,...,n Nf ^. x _ X> ^ x> fa .... x? fa ...) ( 56 ) 

= E n; ( ,7U-:-'« - 1, «*; •••) + Ifo/l-: ^ + 1, a t ; .. 

+ r o ft /(•••; & - 2, /3i; x, f3 2 ; x, /3 m ; ...) 
+ r °fo/(---; ft; x i fa ...,x, /3 m ; ...) 
- n 1 f(...;x,/3 1 ;x,/3 2 ; ...;x,/3 m ; ...), 



and 



-.n-i,ri2,—,n N 



/(...; x, /3 2 ; x, m ; x + l, ft...) (57) 
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= E [ro l a°/(.-; Xi - 1, a i; ...) + r^/(...; x, + 1, a t ; ...) 

i=l,Xj^x,x+l 

+ r oft/(--; ^ ft; ft; -,x, p m ; ...) 

+ r° i V(...;x,/3 2 ;..,x,/3 m ;x + 2,/3 1 ;...) 

-ni/(...;x,/3i;x,/3 2 ; x, ft,; ■■■), 

for neighboring sites occupied by one particle of kind (3\. In (55), (56) and (57) we 
have (3i = 1,2,...,N and = = if /?i 7^ 1, and without loss of generality we 
also choose no collisions of particles 1. Equations (55), (56) and (57) are automatically 
satisfied if /3± 7^ 1. On the other hand, for (3i = 1, while (56) is again automatically 
satisfied, the equations (55) and (57) will impose algebraic constraints for the matrices 
defining the ansatz. By inserting the ansatz (45) with (49) and (50) into equations (55) 
and (57) we obtain, after some algebraic manipulations, the following constraints among 
the matrices: 

4(1) jGfc) AM _ A (M AM 4 (i) ( ro\ 

where fa = 3,4, ...,N and m = 2,3, ... with fa + \ > fa + 1. In order to satisfy (58) the 
for any set {fa, f3 m } we should impose the following commutation relations among 
matrices and A^ a) (a > 3): 

A % A % = A % A % ^ ^ ( 59 ) 

3 3 3 3 

Let us consider now the configurations where we have neighbors sites at positions 
x and x + 1 with multiple particle occupations. For those configurations in which 
21 {I = 1,2,...) particles exchange positions the eigenvalue equation (43) give us the 
following relations: 

e ni>na '-' n "/(...; x, 1; ...; x, 21 - 1; x, ft+i; ...; x, ft*; (60) 
x + 1, 2; x + 1, 2/; x + 1, x + 1, fa m ,; ...) 

= E [r? Q ;/(...; Xi - 1, a,; ...) + r°^/(...; x, + 1, n,: ...)' 

i=l,Xi^!r,:r+l 

+ r oi /(...;x 1, 1; x, 3; x, /3 m ; ...) 
+ r°i/(...;x,2;...;x,2/;x,A+i;-;a:,/3 m ; 

x + 1, 1; x + 1, 2/ - l;x + l,fa l+1 ; ...;x + l,fa m ,; ...) 
- ni/(...;x, 1; ...;x,2Z - l;x, ...;x,(3 m ; 

x + 1, 2; x + 1, 2/; x + 1, x + 1, /? m ,; ...) 

and 

£ ni ' n2 '-' nN /(...; x, 2; x, 2/; x, ft +1 ; ...; x, /3 m ;; (61) 
x + 1, 1; ...;x + 1,2/ - l;x + 1, ft+i; -;a: + l,/3 m ; •••) 

n 

= E [ro l a°/(-; ^ - 1, ",: ...) + r°3/(...; x, + 1, n,: ...) 

i=l,Xi^x,x+l 
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+ 0(...;x + l,3;...;x + lX,;x + 2,l;...) 
+ r o if (...; x, 1; x, 21 — 1; x, x, (3 m r, 

x + 1, 2; x + 1, 2/; x + 1, /3 J+1 ; x + 1, f3 m ; ...) 

— nif(...; x, 2; x, 2/; x, f3 l+1 ; ...; x, /? m/ ; 

x + 1, 1; ...;x + 1,2/ - l;x + 1, A+i! • ••;£ + l,/3 m ; •••)> 

where $+i > 2/ + 1. Furthermore, for those configurations in which 21 + 1 (/ — 0, 1, ...) 
particles exchange positions the eigenvalue equation (43) give us 

e ni > na '-> n "/(...; x, 1; ...; x, 2/ + 1; x, /3 I+2 ; ...; x, /3 m ; (62) 
x + 1, 2; x + 1, 2/; x + 1, x + 1, ...) 

= E 'r;i„7(-:.'-, - 1, n,: ...) + ^/(-l x % + 1, a,; ...)' 

+ Tg if(...]x 1, 1; x, 3; . x, /3 m ; ...) 
+ r?^(...;x,2;...;x,2/;x,A+2;...;a;,/3 m ; 

x + 1, 1; x + 1, 2/ + 1; x + 1, x + 1, ...) 

— ni/(...;x, 1; ...;x,2/ + l;x, A +2 ; ...;x,/3 m ; 

x + 1, 2; x + 1, 21; x + 1, ...; x + 1, p' m ,; ...) 

and 

£ ni '" 2 ' ' JV /(...; x, 2; x, 2/; x, x, /3 m /; (63) 
x + 1, 1; x + 1, 21 + 1; x + 1, $+ 2 ; x + 1, (3 m ; ...) 

= E [rL°/(...; Xi - 1, a,; ...) + VW'Xfi..,. Xl + 1, a t ; .. 

i=l,Xj^x,x+l 

+ T° */(...; x + 1, 3; ...;x + l,/3 m ;x + 2,1;...) 
+ if (■■■', x, 1; x, 2/ + 1; x, x, /3 m/ ; 

x + 1, 2; x + 1, 2/; x + 1, /3 z+2 ; ...; x + 1, f3 m ; ...) 

— n if '(•••! 2, •••) ^, 2/; x, x, /3 m /; 

x + 1, 1; ...;x + 1,2/ + l;x + 1, A+25 x + 1, /3 m ; ...), 

where > 2/ + 2. By inserting the ansatz (45) with (49), (50) and (59) into equations 
(60)-(63) we obtain, after some algebraic manipulations, the following constraints among 
the matrices: 

i(fc(a>+... +fc (M>) .(1) A (2l-1) .(ft+i) 4 (/3 m ) 4 (2) 4 (2i) _ (( , A , 

e A kW ' ' ' ^(^-ij^fcCSf+i) ' ' ' /i ik(/'m) /i fc(2) • • • A k {2t) - yo^) 

i(fc(3) + ... +fc (a-l) )4 (2) 4 (2i) 4 (A+i) 4 (/? m ) ,(1) .(2J-1) 

and 

i(fe (2) + ... +fe (20) 4 (l) 4 (2J+1) 4 (A+i) 4 (/?m) 4 (2) 4 (2«) _ (Rt .. 

6 A kW ' ' ' /1 jfc(2i+l) /i fe (/3 i+ l) ' ' ' /1 fc(,3 m ) /1 fc(2) • • • ^fc(2I) - 

i(fc(3)+...+fc(»+i)) 4 (2) 4 (2«) 4 (ft +1 ) 4 (/3 m ) 4 (1) 4 (2J+1) 

e A km • • ■ /i fc(20 /i fc( ^ +l) • • • /i fc(/3m) /i fc(1) • • • /i fe( 2i+i)- 
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In order to satisfy equations (64) and (65) for all I = 0,1,... and m = 2,3,... with 
A+i > A + 1, the matrices defining the ansatz should satisfy the algebraic relations: 

A a ) AP) - AW A( a ) (R^n-Ll} 

Equations (66) with (49), (52), (53) and (59) completely fix the commutation relations 
among the matrices defining the ansatz: 

EA%=e^A%E, (A%) 2 = 0, (67) 

3 3 \ 3 ' 

Ap) AP) _ G<XP (%.(<*) l.W\AP) 4( Q ) ( n «— 1 O Af\ 

k (a) k^ ~ a P ^ j ' 1 ' kf- ] k {a) ^ ' ^ ~~ ' ' '"' 

where the coupling constants S^(kj a \kf) are given by: 

(1) (1) rjj + rJje^+^-e^- 15 

) = - r lo + i^i c <^)_ c ^' (68) 
^«(*i a) ,*i (a) ) = l (2<«<AT), 

S^ia^S^. = 1 n 7 +T1 = (1< a < iV - 1), 

a+iaV j i i ) oaa+lfAa) h.( a + 1 )\ \ — — ji 

°a a+1 K^l ) ft j J 

<S^(*j a \ ^) = Sfe{kf, kf) = 1 (a = 1, AT - 1, a + 1< /3 < N). 

Finally, all other relations coming from the eigenvalue equation (43) containing 
amplitudes with arbitrary number of particles on neighbors sites are automatically 
satisfied by the ansatz (44) and (45) with (49) and the algebraic relations (67). 
Furthermore, the associativity of the algebra (67) provides a well-defined value 
for any product of matrices and it follows from the fact that the algebra (67) 
is diagonal and the structure constants (68) are c-numbers with the property 
S:S(kfM^ a a (kf\kf) = l(a,0= l,...,iV). 

In order to complete our solution we should fix the spectral parameters {kf} 
(a = 1,...,N). Like in previous section, the algebraic expression in (67) assures that 
any matrix product defining our ansatz can be expressed in terms of a simple matrix 
product A^ 1} ■ ■ ■ A^A^l) ■ ■ ■ A^ 2) ■ ■ ■ A^ } ■ ■ ■ A^ } E L . From the periodic boundary 
condition we obtain: 

aW aW aW a( n ) 



aW ...aW A (l) A (l) ...A {N) E L (69) 

k m k m k w k (1) frW^ ^ Dy J 

K l j-l j K j + 1 Kn N 

-noii/iW tWv^E^i^p-ftfijW ...4(1) A m ... A W F l a (i) 

t> j K l 3-1 J + l Kn N K j 

— Usui's -.h > e q e 3 ^(1)^(1) -^,(1)^.(1) a ,(n)& 



i=i 



1 3-1 3 3+1 n N 
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and 



aW ...A^ A (a) A (a) ...A^ F, L (1C\\ 

K i j-\ j j+i n N 



K l K j-1 K j + 1 K n N 

e 9 e J A t («) A t (i) >)>) 



j K l K j-1 K j + 1 K ™jV 



- e 9 e J V' V» V» V» V 1 ' 



for a = 2, N — 1, and 

K i j-\ j j+i n N 

_ -ilf»i 4 (l) .(JV) j(JV) .(JV) pL 4 (iV) 

~ e ^1.(1) " ' ^- k (N)^- k (N) ■ ^,(N)^ ^,(N) 

1 2 — 1 2 + 1 n iV J 

- p-^r^W 4« . . . A (N) A {N) ■ ■ ■ A (N) E L 

2 1 j — 1 j + 1 n jV 

_ ik f\L-n N ^)Al) AN) AN) AN) AN) pL 

~ G A k W ' ' ' A k W A k (N)^-AN) ■ ■ ■ A ,(N)^ , 

K l 3-1 2 2 + 1 n N 

where in (69), (69) and (71) we used the algebraic relations (67) with (68) and we 



introduced S\\(kf\ k^p) = —1. From (69), (69) and (71) we obtain the Bethe equations 



for our model: 



~ ( ' Mrs! + r;;e< ,+ '." , )- e < > ' ( ' 



e ^ = e -*£ 9= + 4 ( a = 2 , AT - 1), (73) 

e ^ Ar) (L-n JV _ 1 ) _ ^4) 



where equations should be satisfied for all k^ (a = 1, ...,N) with j = 1, ...,n a . On the 
other hand, the momentum of the eigenstate is given by inserting the ansatz (44) and 
(45) into relation (46) and (47): 

n\ n,2 n N n 1 

p = E^ (1) + E^ ) + --- + E^ ) = x C = o.i.-^-i)- (75) 
j=i j=i j=i 

where the spectral parameters satisfy the Bethe equation (72), (73) and (74). 

The Bethe equations (72), (73) and (74) are more complicated than the case of 
previous section since the spectral parameters, and the eigenvalues of our model, depend 
on the densities of impurities. We compare the Bethe equation solution (78) for N = 3 
with the eigenvalues obtained from direct diagonalization of the Hamiltonian (42) with 
L = 5, rii = 2, n2 = 1 and 713 = 1 (see Appendix B). The model displays the full 
spectrum of the ASEP and additional eigenvalues. The spectrum of the ASEP is 
obtained when Ya]=i k^ = for all a = 2, ...,7V. The stationary state belongs to 
this case and has the eigenvalue with the largest real part e ni ' n ' 2 = 0. In order to obtain 
the second largest real part eigenvalue, we can rewritten the Bethe equations in a more 
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convenient way by eliminating the the spectral parameters k^p in equation (72). The 
first excited state is obtained when we have 

V k (a) = — (76) 

^ 3 (L- n N ^){L - n N _ 2 ) ■ ■ ■ (L — n a _0 1 ) 

for all a = 2, N (see Appendix B for a detailed discussion for a small system). Hence, 
by using (76) we can relate the sum over spectral parameters k^p for the first excited 
state with roots of unity: 

n 2 r, N-1 

E^rfena-A*)- 1 . ( 77 ) 

q=l L 1=1 

where Pi — ^ are the densities of particles of kind I = 1, N. Finally, inserting (77) 
into (72) we obtain the following Bethe equation: 

The Bethe equation (78) generalizes [22, 23] and (29) to the case of N — 1 kinds of 
impurities. As in the case of one kind of impurities (N — 2), the phase factor on (78) 
plays a fundamental role in the spectral properties of the model. The eigenvalue with 
the second largest real part determines the relaxation time and the dynamical exponent 
z. This eigenvalue is provided by selecting n\ fugacities from (78) with momentum 
P = ^ in (75). We solve (78) numerically for the totally asymmetric exclusion process 
(N-TASEPI), when rj ° = 1 and T° J = (or ° = and r?J = 1), in the half-filling 
sector ni = L/2 for both two (N = 3) and three (N = 4) kinds of impurities (see 
Figure (4)). As in previous section, the errors displayed are computed from the linear 
regression for the logarithm of the real part of the energy gap versus the logarithm of 
L = 80, 160, 200, 320, 400. For iV = 3 we consider n 2 = L/8 and n 3 = L/8. In this case 

7 3 

we found that the spectral gap scales with L~2 instead of L~2 for the ASEP [6] and 
L - i for the ASEPI [23]. On the other hand, for N = 4 we choose n 2 = L/8, n 3 = L/8 

9 

and n 4 = L/8, and we found that the spectral gap scales with L~a. From these results, 
we can conjecture that for general N = 2,3, ... our model will display a scaling exponent 
z = | + N — 1. This scaling law is a consequence of the dependency of the dynamic of 
impurities of kind N from the dynamic of impurities of kinds a < N, since the impurities 
a only move when particles of kind a — 1 change position with them. As a consequence, 
the time to vanish the fluctuations on the densities of particles of kinds N — 1 acts as 
a time scale for the diffusion of the impurities iV (reflected in the L~5 x L~ N+1 gap for 
the N-ASEPI). 

Eigenstates: The stationary state is the eigenstate associated to the eigenvalue 
£ ni,n 2 — q ? anc [ w j^j 1 mome ntum P = 0. In this case, the Bethe equations (72), (73) 

, (a) 

and (74) have the unique solution e i =1 for all a = 1,...,N and j = l,...,n a , and 
the ^-matrix reduces to the identity S"^°^(kj ai \ki a2 ^) = 1. Consequently, like in the 
previous section, at the stationary state all configurations that satisfy the hard-core 
constraints imposed by the definition of the model can occur with equal probabilities. 
Furthermore, each site is occupied by a particle a with probability p a = n a /L. 
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h- 1 




Ln(L) 



Figure 4. In this figure we display the logarithm of the real part of the energy gap 
versus the logarithm of L for N = 3, 4 in the in the half-filling sector n\ = L/2. The 
dynamical exponents z are 3.50 and 4.01(2), respectively. 



4. Conclusion 



In the present work we formulate an exactly solvable asymmetrical diffusion model of 
N = 1,2,3,... kinds of particles with impurities (N-ASEPI). In this model particles 
of kind 1 can jump to neighboring sites if they are empty and particles of kind a = 
2,3, ...,N (called impurities) only exchange positions with others particles, satisfying a 
well defined dynamics. We solve this model with periodic boundary condition through 
a new matrix product ansatz [24, 21] and we analyze the spectral gap for some special 
cases. Our N-ASEPI model displays the full spectrum of the ASEP [6] plus new 
levels. The first excited state belongs to these new levels and has unusual scaling 
exponents. Although the ASEP belongs to the KPZ universality class, characterized 
by the dynamical exponent z — | [11], we conjecture that our model displays a scaling 
exponent | + N — 1, where N — 1 is the number of kinds of impurities. In order to check 
our conjecture, we solve numerically the Bethe equation with N = 3 and N = 4 for the 
totally asymmetric diffusion and we found that the gap for the first excited state scale 

7 9 

as L~2 and L~2 in these cases. Furthermore, for iV = 2 we generalize the model [23] to 
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include quantum spin chain Hamiltonians and we analyze the Bethe Ansatz equation 
for the symmetric and asymmetric diffusions. A quite interesting problem for the future 
concerns the formulation of the model with open boundary conditions instead periodic 
ones. In this case we expect the critical behavior of the model will display the same 
scaling exponent | + iV — 1 . 
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Appendix A. Eigenvalues and Bethe roots for N = 2, n\ = 2, n 2 = 1, and 

L = 6 

In this appendix we list the full spectrum of the Hamiltonian (3) for the asymmetric 
exclusion model with one kind of impurity for ./V = 2, n\ — 2, n 2 — 1, L — 6, and 
rJ5 = 0.75. We compared the eigenvalues obtained by direct diagonalization of the 
Hamiltonian (3) with those given by the Bethe ansatz solution. For iV = 2, n\ = 2, 
n 2 — 1, L — 6, and Tq? = 0.75 the Bethe equations (32) reduces to 

6 = i¥m 0.75 + 0.25a;y -s 

0.75 + 0.25xy - y ' v ' ; 

where m — 0, 1, 2, 3 and x = e lfc i ) and y = e lfca ) are the fugacities. Equation (A.l) can 
be reduced to a simple polynomial equation by inserting the relation for the momentum 
P = ¥ (31) 

xy = e *¥+'x™ (A.2) 

where I — 0, 1, 2, 3, 4, 5 and we have used (28). In table (Al) we display the full spectrum 
of the Hamiltonian (3) and the associated Bethe roots x and y. The eigenvalue with the 
largest real part is zero, it is provided by choosing m = and P = (/ = 0) in (A.l) and 
(A.2). In this case all spectral parameters are zero and we have e* fc i = e* fc 2 = e ife i = 1. 
It is also important to note that for m = our model reproduces the full spectrum of 
the ASEP. For m — 1,2,3 our model displays additionals energy levels. In special, the 
eigenvalues with the second largest real part, that determine the relaxation time and the 
dynamical exponent z, belong to these new levels. Actually, these eigenvalues are given 
by a complex conjugated pair by choosing m = 1 and P = ^ = ^ and by choosing 
m = L-n 1 - 1 = 3 and P = (L-l)^ = 5f. 

Appendix B. Eigenvalues and Bethe roots for N — 3, n\ — 2, n 2 — 1, n 3 = 1 
and L = 5 



In this appendix we list the full spectrum of the Hamiltonian (42) for the asymmetric 
exclusion model with one kind of impurities for N — 3, rt\ — 2, n 2 — 1, n s — 1, L — 5, 
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Table Al. The eigenvalues for L = 6, n\ = 2 and n 2 = 1 



Sector 


Energies 


Bethe roots 


P = 0, m = 




-1.3819660 
-3.6180340 


i; i 

0.3090170 + 0.9510565* 
-0.8090170 ± 0.5877852i 


P = ^t, m = 


-0.5074856 - 0.2181393i 
-2.4925144 + 0.2181393i 


0.8199801 - 0.0796133i; 0.5024902 + 1.1049419* 
-0.7655328 + 0.7087925*; 0.2122933 - 0.9347130* 


P = 2^f, m = 


-0.8772838 - 0.5974379i 
-1.3562417 - 0.1208056i 
-2.7664745 + 0.7182435* 


0.5755876 + 0.3275559*; -0.0093985 + 1.5099420* 
0.7592145 - 0.0749922*; -0.7637957 + 1.0652413* 
-0.9728436 - 0.4011708*; 0.1255214 - 0.9419611* 


P = 3^, m = 


-1.5 
-1.5 


0.6180340; -1.6180340 
-1.6180340; 0.6180340 


P = 4^, m = 


-0.8772838 + 0.5974379* 
-1.3562417 + 0.1208056* 
-2.7664745 - 0.7182435* 


-0.0093985 - 1.5099420*; 0.5755876 - 0.3275559* 
0.7592145 + 0.0749922*; -0.7637957- 1.0652413* 
0.1255214 + 0.9419611*; -0.9728436 + 0.4011708* 


P = 5^, m = 


-0.5074856 + 0.2181393* 
-2.4925144 - 0.2181393* 


0.8199801 + 0.0796133*; 0.5024902 - 1.1049419* 
0.2122933 + 0.9347130*; -0.7655328 - 0.7087925* 


P = 0, m = 1 


-0.6557872 + 0.2961339* 
-2.3442128 - 0.2961339* 


0.7018547 - 0.0959405*; 0.1911910 - 1.3986613* 
0.4489168 + 0.8140546*; -0.9419626 - 0.5194527* 


P = *t m = 1 
6 ' 


-0.0572751 + 0.2059994* 
-1.3828654 + 0.0130776* 
-3.5598595 - 0.2190769* 


0.8373023 - 0.0855293*; 1.0839923 - 0.4864275* 
0.5293573 + 0.8049461*; 0.0602956 - 1.0362278* 
-0.6226580 + 0.7713741*; -0.9411875 - 0.3629724* 


D O 27T 1 

P — m — l 


-0.4139462 - 0.1150956* 
-2.5860538 + 0.1150956* 


0.6884251 + 0.8429321*; 0.8591734 - 0.3257069* 
-0.5524868 + 0.8570897*; -0.0480098 - 0.9794781* 


P = 3^, m= 1 


-0.4977725 - 0.5315017* 
-1.3835508 - 0.0616583* 
-3.1186767 + 0.5931600* 


0.6345797 + 0.2963134*; 0.6041134 + 1.2937589 
0.7884227 - 0.3472311*; -0.4678531 + 1.0623069* 
-1.0219882 - 0.1463166*; -0.1372745 - 0.9588315* 


P — 4 2-7T — i 


-1.0832361 - 0.2999300* 
-1.9167639 + 0.2999300* 


0.4953389 + 0.1442693*; -1.3406309 + 1.3998737* 
-1.2180653 + 0.0720208*; 0.7326938 - 0.3671649* 


P = 5^, m = 1 


-1.2103358 + 0.1188691* 
-1.4330978 + 0.6646123* 
-2.3565664 - 0.7834814* 


0.5489836 - 0.0641372*; -1.4512935 - 1.0803272* 
0.6425791 - 0.4326605*; -0.5668353 - 1.1597751* 
-0.8735004 + 0.6433986*; 0.3694030 + 0.8445026* 


P = 0, m = 2 


-1.1746854 
-1.9126573 - 0.7734344* 
-1.9126573 + 0.7734344* 


0.4712748; -2.1219040 
0.5609520 + 0.6692321*; -0.7356374 + 0.8776368* 
0.5609520 - 0.6692321*; -0.7356374 - 0.8776368* 


P=^,m = 2 


-0.8494232 + 0.3327914* 
-2.1505768 - 0.3327914* 


0.5692219 - 0.1648546*; -0.4038910 - 1.6383922* 
0.6340571 + 0.6212665*; -1.0851023 - 0.3026347* 


P = 2^f, m = 2 


-0.2257571 + 0.3909776* 
-1.3848207 + 0.0309866* 
-3.3894222 - 0.4219643* 


0.7355889 - 0.2131056*; 0.9417640 - 0.9044865* 
0.6985648 + 0.5996111*; -0.2005790- 1.0675545* 
-0.3934528 + 0.8990944*; -1.0126551 - 0.1129670* 


P = 3^, m = 2 


-0.3819660 
-2.6180340 


0.8090170 + 0.5877853*; 0.8090170 - 0.5877853* 
-0.3090170 + 0.9510565*; -0.3090170 + 0.9510565* 


P = A^f, m = 2 


-0.2257571 - 0.3909776* 
-1.3848207 - 0.0309866* 
-3.3894222 + 0.4219643* 


0.7355889 + 0.2131056*; 0.9417640 + 0.9044865* 
-0.2005790 + 1.0675545*; 0.6985648 - 0.5996111* 
-1.0126551 + 0.1129670*; -0.3934528 - 0.8990944* 


P = 5^p, m = 2 


-0.8494232 - 0.3327914* 
-2.1505768 + 0.3327914* 


0.5692219 + 0.1648546*; -0.4038910 + 1.6383922* 
0.6340571 - 0.6212665*; -1.0851023 + 0.3026347* 


p = o, m = 3 


-0.6557872 - 0.2961339* 
-2.3442128 + 0.2961339* 


0.7018547 + 0.0959405*; 0.1911910 + 1.3986613* 
0.4489168 - 0.8140546*; -0.9419626 + 0.5194527* 


P = ^t, m = 3 


-1.2103358 - 0.1188691* 
-1.4330978 - 0.6646123* 
-2.3565664 + 0.7834814* 


0.5489836 + 0.0641372*; -1.4512935 + 1.0803272* 
0.6425791 + 0.4326605*; -0.5668353 + 1.1597751* 
-0.8735004 - 0.6433986*; 0.3694030 - 0.8445026* 


P = 2^, m = 3 


-1.0832361 + 0.2999300* 
-1.9167639 - 0.2999300* 


0.4953389 - 0.1442693*; -1.3406309 - 1.3998737* 
0.7326938 + 0.3671649*; -1.2180653 - 0.0720208* 


P = 3^, m = 3 


-0.4977725 + 0.5315017* 
-1.3835508 + 0.0616583* 
-3.1186767 - 0.5931600* 


0.6345797 - 0.2963134*; 0.6041134 - 1.2937589 
0.7884227 + 0.3472311*; -0.4678531 - 1.0623069* 
-1.0219882 + 0.1463166*; -0.1372745 + 0.9588315* 


P = 4^, m = 3 


-0.4139462 + 0.1150956* 
-2.5860538 - 0.1150956* 


0.8591734 + 0.3257069*; 0.6884251 - 0.8429321* 
-0.0480098 + 0.9794781*; -0.5524868 - 0.8570897* 


P = 5^, m = 3 


-0.0572751 - 0.2059994* 
-1.3828654 - 0.0130776* 
-3.5598595 + 0.2190769* 


0.8373023 + 0.0855293*; 1.0839923 + 0.4864275* 
0.5293573 - 0.8049461*; 0.0602956 + 1.0362278* 
-0.6226580 - 0.7713741*; -0.9411875 + 0.3629724* 



and TqJ = 1. We compared the eigenvalues obtained by direct diagonalization of the 
Hamiltonian (42) with those given by the Bethe ansatz solution. For N — 3, rt\ — 2, 
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n 2 — 1, ns — 1, L — 5, and I^y = 1 the Bethe equations (72), (73) and (74) reduce to 
x 5 = -e^ k(2) — , (B.l) 

i -y 

e i4fc(3) = 1, (B.3) 

where x = e lk ^ ) and y = e i/c 2 ' are the fugacities. Equation (B.l) can be reduced to a 
simple polynomial equation by inserting the relation for the momentum P = ^ (31) 

xy = e i ™- i % m - i % m ', (B.4) 

where I = 0,1,2,3,4,5, m = 0,1,2,3, w! = 0,1,2, and we have used from (B.2) and 
(B.3) the following relations: 

k^ = -— m+— m', (B.5) 
12 3 

= —m. (B.6) 

In tables (Bl) and (B2) we display the full spectrum of the Hamiltonian (42) and the 
associated Bethe roots x and y. The eigenvalue with the largest real part is zero, it is 
provided by choosing m = 0, m! = and P = (/ = 0) in (B.l), (B.4) and (B.5). In this 
case all spectral parameters are zero and we have e* fc i ' = e lk z ' = e* fc<2) = e* fc<3) = 1. As 
in the case iV = 2, it is important to note that for m = wl = our model reproduces the 
full spectrum of the ASEP. If m or m' are not zero our model displays additionals energy 
levels. In special, the eigenvalues with the second largest real part, that determine the 
relaxation time and the dynamical exponent z, belong to these new levels. Actually, 
these eigenvalues are given by a complex conjugated pair by choosing m — 1, m' — 
and P = y- = \ and by choosing m = L-n 2 -l = 3, m' — 1 and P = (L - 1) ^ = 4^. 
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Table Bl. The eigenvalues for L = 5, n\ — 2, n 2 — 1 and n 3 = 1 



Sector 


Energies 


Bcthe roots 


P = 0, m = 0, rri = 




-2 


i; i 

i; —i 


P = ^L, m = 0, m' = 


-0.7102901 + 0.3268955* 
-3.0987269 - 0.9146808* 


0.7044758 + 0.2260418i; 0.0049626 - 1.3516126i 
-0.1522713 + 0.9227363i; -1.0571671 - 0.1604367i 


P = 2^f, to = 0, to' = 


-1.2113723 + 1.4418449* 
-1.4796107- 0.4907884i 


-0.4434074 - 1.1954323i; 0.6528893 - 0.4345885i 
0.6498509 + 0.3220402i; -1.3593328 - 0.2308612* 


P = 3^, m = 0, m' = 


-1.2113723 - 1.4418449* 
-1.4796107 + 0.4907884* 


-0.4434074 + 1.1954323i; 0.6528893 + 0.4345885i 
0.6498509 - 0.3220402i; -1.3593328 + 0.2308612i 


P = 4^, to = 0, to' = 


-0.7102901 - 0.3268955* 
-3.0987269 + 0.9146808* 


0.7044758 - 0.2260418i; 0.0049626 + 1.3516126i 
-0.1522713 - 0.9227363i; -1.0571671 + 0.1604367i 


P = 0, m = 0, to' = 1 


-0.9146286 + 0.4740955i 
-2.5853714- 1.3401209* 


0.6805844 + 0.0435103i; -0.8126913 - 1.2205172* 
0.2122293 + 0.8634034i; -1.0801224 + 0.3136036* 


P=?£,m, = 0, m' = l 


-0.1515123 + 0.6102038* 
-1.9349423 - 0.2034672* 


0.7626883 - 0.0415504*; 0.9276613 - 0.9238376i 
0.3562835 + 0.8258711*; -0.4639570 - 1.0103645* 


P = 2^f, to = 0, to' = 1 


-0.5999763 - 0.1164667* 
-3.3781712 + 0.3243784* 


0.5948634 + 0.9609114*; 0.7314932 - 0.4978679* 
-0.8574914 + 0.5560440*; -0.5334672 - 0.8202624* 


P = 3^, to = 0, to' = 1 


-0.5812538 - 1.117557H 
-1.7496157 + 0.3744123i 


0.7168460 + 0.2139355*; 0.2462898 + 1.3138552i 
-0.9667610 + 0.7800852i; 0.5682275 - 0.5702092* 


P = 4^E, TO = 0, TO' = 1 


-1.1838852 - 0.5304177i 
-1.9206433 + 1.5249396i 


0.6689977 + 0.1247276*; -1.3569984 + 0.5637791* 
-0.8877912 - 0.7977974*; 0.4931159 - 0.6773195* 


P = 0, to = 0, m' = 2 


-0.9146286 - 0.4740955i 
-2.5853714 + 1.3401209i 


0.6805844 - 0.0435103*; -0.8126913 + 1.2205172* 
-1.0801224 - 0.3136036*; 0.2122293 - 0.8634033* 


P = ^t, m = 0, m' = 2 


-1.1838852 + 0.5304177i 
-1.9206432 - 1.5249396i 


0.6689977 - 0.1247276*; -1.3569984 - 0.5637791* 
0.4931159 + 0.6773194*; -0.8877913 + 0.7977974* 


P = 2^, to = 0, m' = 2 


-0.5812538 + 1.1175571i 
-1.7496155 - 0.3744122i 


0.7168460 - 0.2139355*; 0.2462898 - 1.3138552* 
0.5682275 + 0.5702092*; -0.9667610 - 0.7800853* 


P = 3^, to = 0, m' = 2 


-0.5999764 + 0.1164667i 
-3.3781713 - 0.3243784i 


0.7314931 + 0.4978680*; 0.5948633 - 0.9609114* 
-0.5334673 + 0.8202623*; -0.8574914 - 0.5560439* 


P = 4^, m = 0, to' = 2 


-0.1515123 - 0.6102038i 
-1.9349423 + 0.2034672i 


0.7626883 + 0.0415504*; 0.9276613 + 0.9238376* 
-0.4639570 + 1.0103645*; 0.3562835 - 0.8258711* 


P = 0, m = 1, to' = 


-0.6730048 + 0.2787791i 
-3.1930207- 0.7787791i 


0.7126887 + 0.2811734*; 0.1922387- 1.2909955* 
-0.2488952 + 0.9128580*; -1.0220576 - 0.2690613* 


P = m = 1, to' = 


-0.0095876 - 0.1567928i 
-1.9958904 + 0.0522644i 


0.8542110 - 0.0537788*; 1.1253051 + 0.3142423* 
-0.1080056 + 1.0208495*; 0.1011590 - 0.9688728* 


P = 2^e, to = 1, to' = 


-0.7535156 - 0.3709970i 
-2.9896292 + 1.0401277i 


0.6970108 - 0.1760220*; -0.1977531 + 1.3768984* 
-1.0811538 + 0.0468434*; -0.0567203 - 0.9223284* 


P = 3^, to = 1, to' = 


-1.3863511 - 1.4870626i 
-1.4057371 + 0.5089149i 


0.6242368 + 0.4958318*; -0.5799899 + 1.1122610* 
0.6576010 - 0.2680371*; -1.4074815 + 0.0448287* 


p = 4^f, to = 1, to' = 

o 


-1.0413131 + 1.3814181i 
-1.5519502 - 0.4678727i 


0.6752325 - 0.3750635*; -0.2901255 - 1.2617289* 
0.6383338 + 0.3796151*; -1.2858348 - 0.3995129* 


P = 0, to = 1, m' = 1 


-1.2570659 - 0.5290855i 
-1.7429341 + 1.5290855i 


0.6663053 + 0.1698506*; -1.4092394 + 0.3592349* 
-0.8018074 - 0.9104995*; 0.5447415 - 0.6185860* 


P = ^f, m = 1, to' = 1 


-0.8561485 + 0.4450459i 
-2.7316367 - 1.2540628i 


0.6850761 + 0.0857982*; -0.6152818 - 1.3111923* 
0.1266103 + 0.8921177*; -1.0932067 + 0.1912369* 


P = 2^t, TO = 1, to' = 1 


-0.0857230 + 0.4635491i 
-1.9632204 - 0.1545321i 


0.7818977 - 0.0047923*; 1.0392523 - 0.7453724* 
0.2799188 + 0.8805384*; -0.3409951 - 1.0271761* 


P = 3^, TO = 1, to' = 1 


-0.6175870 - 0.1729136i 
-3.3334697 + 0.4819306i 


0.4918745 + 1.0873199*; 0.7281574 - 0.4146480* 
-0.9214012 + 0.4677083*; -0.4406702 - 0.8616120* 


P = 4^, TO = 1, to' = 1 


-0.7242529 - 1.2183559i 
-1.6879620 + 0.4093389i 


0.7058518 + 0.2643779*; 0.0586461 + 1.3254225* 
-1.0844991 + 0.6747096*; 0.5987696 - 0.5044365* 


P = 0, to = 1, m' = 2 


-0.2349507 - 0.7502975i 
-1.8990239 + 0.2502975i 


0.7483302 + 0.0808009*; 0.7839712 + 1.0726279* 
-0.5896583 + 0.9786229*; 0.5447415 - 0.6185860* 


P = m = 1, to' = 2 


-0.9772050 - 0.4973541i 
-2.4295315 + 1.4108995i 


0.6769329 - 0.0020008*; -0.9917103 + 1.0948804* 
-1.0535823 - 0.4366956*; 0.2924923 - 0.8265846* 


P = 2^, to = 1, to' = 2 


-1.1124249 + 0.5256078i 
-2.0954867 - 1.5037553i 


0.6714067 - 0.0814666*; -1.2684629 - 0.7597094* 
0.4335917 + 0.7324188*; -0.9579927 + 0.6801679* 


P = 3^, to = 1, to' = 2 


-0.4510354 + 1.0052159i 
-1.8058199 - 0.3360853i 


0.7267522 - 0.1667177*; 0.4348678 - 1.2686880* 
0.5291823 + 0.6366906*; -0.8431292 - 0.8649373* 


P = 4^, to = 1, to' = 2 


-0.5893425 + 0.0585923i 
-3.4051795 - 0.1631208i 


0.7266229 + 0.5951587*; 0.6653905 - 0.8311391* 
-0.6226556 + 0.7686536*; -0.7857321 - 0.6360565* 



[13] C. Boldrighini, G. Cosimi, S. Frigio, and G. Nunes, J. Stat. Phys. 55, 611 (1989); P. A. Ferrari, 
C. Kipnis and E. Saada, Ann. Prob. 19, 226 (1991), P. A. Ferrari, Prob. Theory Relat. Fields 
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Table B2. The eigenvalues for L = 5, n\ — 2, n 2 — 1 and n 3 = 1 



Sector 


Energies 


Bcthe roots 


P = 0, m = 2, to' = 


-0.8782556 + 1.3066224i 
-1.6217444 - 0.4405970i 


0.6924404 - 0.3181480i; -0.1217444 - 1.3066224i 
0.6217444 + 0.4405970i; -1.1924404 - 0.5478775* 


P = ^L, m = 2, m' = 


-0.6420067 + 0.2272404i 
-3.2715388 - 0.6339770* 


0.7210291 + 0.3432755i; 0.3565184 - 1.2004077i 
-0.3454212 + 0.8925020i; -0.9765410 - 0.3717780i 


P = 2^f, m = 2,m'=0 


-0.0382564- 0.3118707i 
-1.9835959 + 0.1039590i 


0.8094973 - 0.0284654i; 1.1094939 + 0.5414704i 
-0.2219898 + 1.0302847i; 0.1946917 - 0.9286414i 


P = 3^, m = 2, to' = 


-0.8022829 - 0.4105415i 
-2.8668478 + 1.1536863i 


0.6905308 - 0.1296573i; -0.4074337 + 1.3637265i 
-1.0933657 - 0.0707271i; 0.0366096 - 0.9119650i 


P = 4^, m = 2, m' = 


-1.3312716 + 0.5218134i 
-1.5641999 - 1.5163354i 


0.6627602 - 0.2174407*; -1.4253665 - 0.1539334i 
0.5883753 + 0.5576163i; -0.6993882 + 1.0161915* 


P = 0, m = 2, to' = 1 


-0.8782556 - 1.3066224i 
-1.6217444 + 0.4405970i 


0.6924404 + 0.3181480*; -0.1217444 + 1.3066224* 
0.6217444 - 0.4405970*; -1.1924404 + 0.5478775* 


P=^,m = 2, to' = 1 


-1.3312716 - 0.5218134i 
-1.5641999 + 1.5163354i 


0.6627602 + 0.2174407*; -1.4253665 + 0.1539334* 
-0.6993882 - 1.0161915i; 0.5883753 - 0.5576162* 


P = 23r, m = 2, m' = 1 


-0.8022829 + 0.4105415i 
-2.8668476 - 1.1536863i 


0.6905308 + 0.1296573*; -0.4074337 - 1.3637265* 
0.0366097 + 0.9119649*; -1.0933658 + 0.0707270* 


P = 3^, m = 2, m' = 1 


-0.0382564 + 0.3118707i 
-1.9835961 - 0.1039590i 


0.8094973 + 0.0284654*; 1.1094939 - 0.5414704i 
0.1946916 + 0.9286414*; -0.2219899 - 1.0302847* 


P = 4^e, m = 2, m' = 1 


-0.642007 - 0.2272403i 
-3.2715388 + 0.6339770i 


0.3565184 + 1.2004077*; 0.7210291 - 0.3432755* 
-0.9765410 + 0.3717780*; -0.3454213 - 0.8925020* 


P = 0, m = 2, to' = 2 


-0.5857864 
-3.4142136 


0.7071068 + 0.7071068*; 0.7071068 - 0.7071068i 
-0.7071068 + 0.7071068*; -0.7071068 - 0.7071068* 


P = m = 2, to' = 2 


-0.3351580 - 0.8824028i 
-1.8558250 + 0.2946176i 


0.7368286 + 0.1224476*; 0.6168508 + 1.1882338* 
-0.7166508 + 0.9306243i; 0.4810054 - 0.7024638* 


P = 2^, m = 2, m' = 2 


-1.0433294 - 0.5145741i 
-2.2656875 + 1.4656306i 


0.6739498 + 0.0394376*; -1.1454536 + 0.9391785* 
-1.0130569 - 0.5593159i; 0.3665267 - 0.7825715* 


P = 3^, m = 2, m' = 2 


-1.0433294 + 0.5145741i 
-2.2656875 - 1.4656306i 


0.6739498 - 0.0394376*; -1.1454536 - 0.9391785* 
-1.0130569 + 0.5593159i; 0.3665267 + 0.7825715* 


P = 4^, m = 2, m' = 2 


-0.3351580 + 0.8824028i 
-1.8558250 - 0.2946176i 


0.7368286 - 0.1224476*; 0.6168508 - 1.1882338* 
0.4810054 + 0.7024638*; -0.7166509 - 0.9306242* 


P = 0, m = 3, to' = 


-1.2570659 + 0.5290855i 
-1.7429341 - 1.5290855i 


0.6663053 - 0.1698506*; -1.4092394 - 0.3592349* 
0.5447415 + 0.6185860*; -0.8018074 + 0.9104995* 


P = m = 3, to' = 


-0.7242529 + 1.2183560i 
-1.6879619 - 0.4093389i 


0.7058518 - 0.2643779*; 0.0586461 - 1.3254224* 
0.5987696 + 0.5044365*; -1.0844991 - 0.6747097* 


P = 2^E, m = 3, to' = 


-0.6175868 + 0.1729136i 
-3.3334696 - 0.4819305i 


0.7281573 + 0.4146480*'; 0.4918745 - 1.0873199i 
-0.4406701 + 0.8616120*'; -0.9214011 - 0.4677083* 


P = 3nr, to = 3, to' = 


-0.0857230 - 0.4635491i 
-1.9632206 + 0.1545321i 


0.7818977 + 0.0047922*'; 1.0392523 + 0.7453723* 
-0.3409952 + 1.0271760*; 0.2799187 - 0.8805385* 


P = 4^: m = 3, m' = 

o 


-0.8561485 - 0.4450459i 
-2.7316368 + 1.2540628i 


0.6850761 - 0.0857982*; -0.6152818 + 1.3111923* 
-1.0932067 - 0.1912369i; 0.1266102 - 0.8921177* 


P = 0, m = 3, to' = 1 


-0.6730048 - 0.2787791i 
-3.1930206 + 0.7787790i 


0.7126887 - 0.2811734*'; 0.1922387 + 1.2909955* 
-1.0220576 + 0.2690613*'; -0.2488952 - 0.9128580* 


p=^L im = 3 im ' = l 


1.0413131 - 1.3814181i 
-1.5519502 + 0.4678727i 


0.6752325 + 0.3750635*; -0.2901255 + 1.2617289* 
-1.2858348 + 0.3995130i; 0.6383338 - 0.3796150* 


P = 2^t, to = 3, to' = 1 


-1.4057371 - 0.5089149i 
-1.3863511 + 1.4870625i 


0.6576010 + 0.2680371*; -1.4074815 - 0.0448287* 
-0.5799900 - 1.1122609i; 0.6242368 - 0.4958318* 


P = 3^, to = 3, to' = 1 


-0.7535156 + 0.3709970i 
-2.9896293 - 1.0401277i 


0.6970108 + 0.1760220*; -0.1977531 - 1.3768984* 
-0.0567204 + 0.9223285*; -1.0811538 - 0.0468433* 


P = 4^e, TO = 3, to' = 1 


-0.0095876 + 0.1567928i 
-1.9958905 - 0.0522644i 


0.8542110 + 0.0537788*; 1.1253051 - 0.3142423* 
0.1011589 + 0.9688728*; -0.1080056 - 1.0208495* 


P = 0, to = 3, to' = 2 


-0.2349507 + 0.7502975i 
-1.8990239 - 0.2502976i 


0.7483302 - 0.0808009*; 0.7839712 - 1.0726279i 
0.4233823 + 0.7660263*; -0.5896583 - 0.9786229* 


p = m = 3, m ' = 2 


-0.5893423 - 0.0585923i 
-3.4051794 + 0.1631208i 


0.6653905 + 0.8311391*'; 0.7266229 - 0.5951587i 
-0.7857321 + 0.6360565*; -0.6226555 - 0.7686536* 


P = 2?f, to = 3, to' = 2 


-0.4510354 - 1.0052159i 
-1.8058198 + 0.3360852i 


0.7267522 + 0.1667177*'; 0.4348678 + 1.2686880* 
0.5291824 - 0.6366905*; -0.8431291 + 0.8649374* 


P = 3^, m = 3, m' = 2 


-1.1124249 - 0.5256078i 
-2.0954868 + 1.5037554i 


0.6714067 + 0.0814665*; -1.2684629 + 0.7597094* 
-0.9579926 - 0.6801680i; 0.4335916 - 0.7324189* 


P = 4^, to = 3, to' = 2 


-0.9772050 + 0.4973541i 
-2.4295316 - 1.4108996i 


0.6769329 + 0.0020008*; -0.9917102 - 1.0948804* 
0.2924923 + 0.8265847*; -1.0535823 + 0.4366957* 
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